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An  e-.'l'er  explanation  of  sunrise  and  sunset  fading  on  long  VLF  paths  is  confirmed  by  examining 
^  data  taken  ov^r  a  wider  range  of  frequency  and  directions  of  propagation.  It  is  found,  however,  that 

J*r»  J  there  is  an  appreciable  dependence  of  the  fading  period  on  the  magnetic  direction  of  propagation.  It 

is  shown  that  this  dependence  is  qualitatively  in  agreement  with  what  would  be  expected  from  a  sharply 
bounded  ionospl.ere  and  a  transverse  magnetic  field. 


1.  Introduction 


In  an  earlier  paper  [Croinbie,  1964]  (subsequently 
referred  to  as  paper  1),  some  of  the  characteristic 
features  of  the  remarkable  periodic  fading  of  VLF 
signals  observed  over  two  long  paths  were  described. 
The  directions  of  these  two  paths  (Canal  Zone  to 
Germany,  18  kc/s,  Hawaii  to  Boulder,  Colo.,  19.8  kc/s) 
were  both  mainly  to  the  east.  On  the  basis  of  the 
data  obtained  from  t*“  tm  it  was  concluded  that,  at 
sunrise,  the  major  features  of  the  observt  1  fading 
could  be  interpreted  as  being  essentially  due  to  inter¬ 
ference  between  the  first-order  modes  propagating  in 
the  daytime  earth-ionosphere  waveguide.  It  was  sug¬ 
gested  that  these  two  modes  were  the  result  of  con¬ 
version  '  of  the  first-  and  second-order  modes 
propagating  in  the  nighttime  portion  of  the  waveguide. 
Since  the  latter  modes  had  different  phase  velocities, 
the  relative  phase,  and  hence  the  combined  amplitude 
and  phase  of  the  two  resulting  first-order  modes, 
depended  on  the  distance  of  the  converting  region 
(possibly  the  sunrise  line)  from  the  source  of  the 
first-  and  second-order  modes,  and  not  on  the  distance 
of  the  receiver  from  the  sunrise  line.  It  seemed 
probable  that  the  transmitting  antenna  v;as  the  source 
of  the  two  initial  modes. 

At  sunset,  however,  for  the  particular  geography 
involved,  it  was  necessary  that  the  (assumed  single) 
mode  incident  in  the  sunset  boundary  be  converted 
at  the  boundary  into  two  modes  of  the  first  and  second 
order.  These  subsequently  interfere  in  the  nighttime 
region  to  the  east  of  the  sunset  line,  causing  the  sunset 
fading.  In  this  case  the  relative  amplitude  of  the  two 
modes  depends  on  the  distance  of  the  receiver  from 
the  boundary. 


1  This  work  wis  supportrd  in  p»rt  by  thr  Aihanonl  Brsrtrch  Project*  Agency.  Wish- 
inglon,  D.C.,  under  Order  No.  i83,  «nd  by  the  National  Science  Foundation  under  GP-.KW2. 


The  feature  common  to  both  the  sunrise  and  sunset 
effects  was  the  presence  of  two  significant  modes  in  the 
nighttime  portion  of  (he  guide,  and  that  the  second-order 
r.iodes  v  is  weaker  than  the  first,  at  the  frequencies 
used.  F  thennore,  in  both  cases  the  characteristics 
of  the  fading  were  determined  by  the  interference  of 
first-  and  second-order  modes  and  the  movement  of 
the  sunrise  or  sunset  line  (the  “terminator”)  along  the 
path.  As  the  terminator  moved  from  one  position 
where  the  two  modes  were  in  antiphase  to  the  next, 
a  distance  D,  the  received  signal  passed  through  one 
complete  cycle  of  the  interference  pattern.  The 
distance  D  is  thus  given  by 


where  At  and  \>  were  the  guide  wavelengths  of  the 
first-  and  second-order  modes  in  the  nighttime  portion 
of  the  earth-ionosphere  waveguide.  By  using  a  simple 
sharply  bounded  isotropic  model  of  the  ionosphere  to 
determine  At  and  A2  it  was  found  that  the  appropriate 
ionosphere  height  was  about  85  km,  confirming  that 
the  mode  interference  was  occurring  in  the  nighttime 
portion  of  the  path.  Other  earth-ionosphere  waveguide 
models  would  of  course  have  given  slightly  different 
heights. 

Although  the  data  used  in  the  above  work  yielded 
a  consistent  picture,  they  pertained  to  a  rather  limited 
situation  since,  as  mentioned  above,  the  two  frequen¬ 
cies  were  close  together  and  the  paths  were  in  similar 
directions.  Nevertheless,  for  each  path  about  1  year 
of  data  was  available  and  used. 

Since  the  earlier  paper  was  written,  the  author  has 
had  an  opportunity  to  exam:ne  more  data  of  limited 
duration  but  covering  a  much  wide"  range  of  frequen¬ 
cies  and  directions  of  propagation.  This  paper  will 
be  concerned  therefore  in  demonstrating  the  extent 
to  which  the  new  data  confirm  the  earlier  conclusions, 
and  in  accounting  for  the  difference. 
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vel'Kiity,  mode  conversion. 
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2.  Experimental  Data 

Figure  l  shows  much  of  the  new  data.  For  each 
path  examined,  the  distance  l)  was  determined  from 
geometry  [Urary,  1%B]  and  I ht  tffw*  rv^ti  Time  between 
signal  minima  (see  paper  1).  The  value  of  1)  was  then 
pbilicj  fUe  ipspfwpfjute  The  jmirM* 

NPM-Bo  at  19.8  ke/s  and  NBA-Fr  | Brady  et  al„  1963| 
at  18  kc/s  arc  the  mean  of  the  observations  given  in 
pa|R  i  1.  flic  ^lata  at  ll)  kl./s  ea  till  NllA-tbi  and 
NBA-Maui  paths  have  been  taken  from  previously 
(uiltiblh'd  (IvitiCh  i  i  all  I  'HH*.  U]  Tlu-  2i~Jti4* 

points  have  been  taken  from  1  month's  observations 
on  NBA  at  Franklurt  IFit/.cnbcrger,  private  communi¬ 
cation!,  Maui  {Katahara,  private  eommunieationj,  and 
Boulder  (Steele.  private  communication].  The 
WWVL-Maui  data  are  from  Morgan  and  Blair  [1965], 


FREQUENCY,  kc/s 


Km. IRK  1,  Effect  of  frequency  and  direction  of  propagation  on  the 
distance,  I),  between  minima  at  sunrise, 

The  small  arrows  ut  v,n  li  point  «ltm*  qiutliltitiw'h  the  rlli'rlitt*  fitrerlinti  of 
with  ri,«|»,il  (•>  thr  h«»ri/on'al  o.injmiirn!  of  tin*  earth's  mngin-tii  fu  hl  (VlagnHir  North 
i-  at  »*!♦’  top  «  f  the  |i.iue.i  The  ilii'hctl  line*  slum  the  \ahtes  of  if  which  would  he  r\|»»  et»  *( 
il  the  ionosphere  ronduetiwtv  gradient  w,i«  O.o/kui  and  for  n  lerenrt'  heights  of  <#)  km 
U't»r\e  1 1  and  70  km  (tun*1  2k 

the  NPM-Tnkvo  point  from  Nakajima  et  al.  |I963|. 
while  the  NBA-New  Zr aland  point  is  from  the  paper 
by  Burtt  (1963).  Observations  at  Stanford  Univer¬ 


sity  on  NSS  (21.4  ke/s)  also  are  in  conformity  with  the 
data  shown  in  figure  1,  although  they  are  not  shown. 

'file  small  arrows  attached  to  each  point  in  figure  I 
are  intended  to  show  cpiulitativcly  the  average  hearing 

id  tilt  path  With  11  S|lci  t  i<  tilt,  ill  1  It  1 1.  *1 1  id  tht  fit  If  I - 

/.on let  component  of  the  earth's  magnetic  field.  This 
ftVr'Tapf'  Pt'Amig  T7  '.tt  tglrted  ir  leHrrf  -Ol' the  fnVt-Cl  (vti 
tude  portions  of  the  paths  where  the  horizontal  com¬ 
ponent  of  the  field  is  greatest. 

Figure  1  shows  Hvu  iiutio  t Hi  t  Is.  that  id  fri-quelicy 
on  /)  and  that  of  the  direction  of  propagation  on  I). 
‘flies*  *,iU  hr.**  U*  *l*seH**ed. 


3.  Effect  of  Frequency 

The  data  used  in  paper  1  were  obtained  for  frequen¬ 
cies  of  18  and  19.8  ke/s.  The  lange  of  frequencies 
shown  in  figure  I  is  from  18  tit  24  kc/s.  | lit*  increase 
resulting  mainly  from  the  change  in  frequency  of  NBA 
from  18  to  24  kc/s.  This  change  produced  relatively 
little  effect,  the  observations  it  Frankfurt  and  Boulder 
showing  a  small  increase  in  U  with  the  increase  in 
frequency,  while  at  Maui  there  was  essentially  no 
change.  It  should  he  noted  that  although  l)  was 
relatively  unchanged  by  the  change  in  frequency,  the 
shape  iff  th»'  sunrise  phase  curves  was  considerably 
altered.  This  is  due  to  the  much  Micreased  amplitude 
of  the  second-order  mode  which  ai  times  exceeds 
that  of  the  first-order  mode.  As  a  result  (see  paper  1) 
the  direction  of  the  rapid  phase  change  at  the  signal 
minima  is  reversed.  If  we  exclude  the  four  points 
for  which  D  >  2400  km,  shown  in  figure  1.  w  hich  will 
be  discussed  later,  it  is  clear  that  the  remainder  are 
in  general  agreement  with  the  behavior  shown  by 
the  dashed  lines.  These  lines  show  the  values  of 
l)  computed  from  phase  velocities  given  bv  Wait  and 
Spies  (1964]  for  isotropic  ionosphere  models  which 
can  be  represented  in  the  form 

fa),  =cu,(/i)  exp  F/3(/i  —  2)], 

where  o>,  =  ofil i>  at  the  height  ;  and  It  is  the  height  at 
which  —  2.5  X  1 CF’:  v  and  are  the  collision  fre¬ 
quency  and  angular  plasma  frequency  at  the  height 
z.  The  lines  in  figure  1  are  for  /i~80  an  1  90  km  and 
/3  —  0.5.  As  noted  above,  they  are  in  general  agree¬ 
ment  w  ith  the  experimental  points  lie  low  /4  =  2400  km. 

4.  Effect  of  Direction  of  Propagation 

Figure  1  apparently  shows  a  very  pronounced  de¬ 
pendence  of  I)  on  the  direction  of  propagation  with 
respect  to  the  direction  of  the  earth’s  magnetic  ft*  Id. 
II  we  divide  the  observations  into  three  groups  of 
direction,  namely,  the  quadrants  centered  on  the  east 
(NPM-Bouldcr,  NBA-Fraukfirrt),  on  the  north  or  south 
(NBA  Bouldcr).  and  on  the  west  (WWVL-Maui,  NBA 
to  Maui,  NPM-Tokyo  and  NBA-New  Zealand),  and 
take  the  mean  value  of  D  for  each  group  irrespective 
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til  frequency,  tlit*  results  ean  he  represented  as  in 
talde  1.  Despite  ihe  conHieting  WVVVL-Maui  res  tilt, 
it  ean  he  seen  that  tliere  is  a  strong  dependence  on  the 
direction  of  propagation,  '  his  out*  discordant  value, 
W  \Y\  1,  (in  Colorado)  to  Maui,  shown  in  figure  1,  is  for 
a  path  which  has  a  much  greater  proportion  (about 
one  third)  of  land  along  it  than  any  of  the  other  paths. 
This  does  not  explain  the  anomalous  value,  however. 
Although  there  is  little  difference  in  the  values  of  I) 
foi  propagation  to  the  east  and  to  the  north  or  south, 
the  valu»“  of  I)  for  propagation  to  the  west  is  appreci¬ 
ably  greater. 

5.  Source  of  Second-Order  Mode 

Figures  2  and  3  show  the  average  rates  of  change 
of  phase  of  the  18-ke/s  signals  from  NBA  as  observed 
at  Maui  (Bradv  et  al.,  lO-Kih]  at  sunrise  and  sunset. 
The  geographical  relation  of  the  transmitter.  rcc<  iver. 
and  terminator  is  also  shown.  At  sunrise  on  this  path 
(which  occurs  between  10  and  If)  LT  for  tin*  months 
shown  m  Itg.  2).  dii  ii<  pin  ol  fading  lUCiraSe”  a.-'  tin 
terminator  approaches  the  receiver.  It  is  evident  that 
any  second-order  mode  excited  by  the  transmitter  is 
being  propagated  towards  the  receiver  over  increasing 
lengths  of  illuminated  path  in  which  the  second  mode 


Taiu.k  1 

thrcrli'Wi  Mean  / )  iinn 

To  the  eaM  ,  1W0 

I'n  the  north  or  south.. 

TolhivuM  im) 


is  highly  attenuated.  Thus  the  increased  depth  of 
fading  observed  as  the  terminator  approaches  the 
receiver  is  a  result  of  the  second-order  mode  generated 
in  tin-  vicinity  of  the  terminator.  Figure  3  lor  the  same 
path  at  sunset  (which  commences  at  the  'ransmitter 
between  00  IT  and  02  l  T)  shows  the  effect  of  the 
second-order  mode  generated  at  the  transmitter. 
The  depth  of  fading  is  greate  t  when  the  terminator 
is  closest  to  the  transmitter  and  thus  represents  the 
effect  of  the  second-older  mode,  excited  by  tile  trans¬ 
mitter.  which  is  converted  into  a  first-ord'*r  nimh*  al 
the  terminator. 

These  observations  clarify  an  issue  which  could  not 
Iw*  ldt-'ek4ed  ifi  pafer  f,  ikat  *.  w-hHiler  nor  il  m*s 
necessary  for  the  second-order  mode  to  be  excited 
only  by  the  transmitter.  It  is  clear  from  the  observa¬ 
tions  at  Maui  that  this  is  not  necessary,  and  that  the 
second-order  mode  which  is  present  in  the  nighttime 


TIME,  UT 


TIME,  UT 


R  T 

SUNRISE 


NBA  -MAUt 
18  kc/s 

Kiel  RE  2.  Holt'  nf  r han fir  ofpha.se  observed  at  Maui  on  the  IR-kc/s 
si/inals  from  \RA  al  sunrise. 


R  T 

SUNSET 


NBA -MAUt 
18  kc/s 

Kiel  RE  3.  Hale  of  eh  an  fie  of  phase  observed  al  Maui  on  llie  I  R-kct 
si  fin  a  Is  from  SHS  al  .sunset. 
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portion  of  the  path  can  also  be  produced  at  or  near  the 
terminator. 


6.  Discussion:  Effect  of  Direction  of 
Propagation 

The  most  striking  feature  of  figure  1  is  the  depend¬ 
ence  of  D  on  the  direction  of  propagation.  This 
suggests  rathei  strongly  thai  il  may  be  due  to  a 
nonreciprocal  effect  on  the  difference  in ‘phase  veloc¬ 
ities  of  the  first-  and  second-order  modes. 

Let  us  consider  an  idealized  earth-ionosphere  wave¬ 
guide  in  which  the  earth  and  ionosphere  are  flat  and 
perfectly  reflecting,  but  the  ionosphere  changes  the 
phase  of  the  reflected  wave  by  6  radians.  The  guide 
wavelength,  of  the  nth  mode  of  a  wave  whose 
free-space  wavelength  is  Kt  is 


K  =  Vl-C?„  (2) 

where 

Cn—fn  —  Gj2n)  (Xol2h),  (3) 

and  h  is  the  height  of  the  uniform  ionosphere  which  is 
assumed  to  have  a  sharp  lower  boundary.  The 

distance  D=  -  ^ ,  is  then  given  by 

—  At  K->  A  i 

D=ma-Q)(  W2), 


since 

1  /  Vl”-  cs  =  1  +  \  c% 

Thus 

D  =  (8/i2/\o)  /  .3  —  OItt).  (4) 

It  is  worth  noting  that  when  6  =  tt,  a  reasonable  ideali¬ 
zation  for  a  sharply  bounded  isotropic  ionosphere, 

a  result  obtained  earlier  [Crombie,  1962].  Equation 
(4)  shows  that  if  D  is  to  he  increased  by  variations  in 
ft .  then  6  must  increase.  To  the  extent  to  which  this 
idealized  model  represents  the  actual  case,  the  de¬ 
pendence  of  D  on  the  direction  of  propagation  must 
result  from  a  value  of  ft ,  the  phase  shift  on  reflection, 
which  is  greater  for  east-to-west  propagation  than  for 
propagation  in  the  opposite  direction.  This  is  just  the 
behavior  found  by  Barber  and  Crombie  [1959]  for  a 
sharply  bounded  ionosphere,  and  thus  qualitatively 
accounts  for  the  observations. 


Equation  (3)  is  a  rather  crude  approximation  since 
losses  in  the  earth  and  ionosphere  have  been  neglected. 
As  a  result  the  guide  wavelength  as  calculated  from 
(2).  and  the  corresponding  phase  velocity  which  is 
proportional  to  \H,  might  be  in  error.  More  exact 
calculations  for  a  sharply  bounded  ionosphere  have 
shown  [Wait,  1962],  however,  that  the  phase  velocity 
of  the  first-order  mode  is  less  for  east-to-west  propaga¬ 
tion  than  for  the  opposite  direction.  This  supports 
the  conclusion  reached  above  since  (2)  and  (3)  show 
that  Xn,  and  hence  the  phase  velocity,  is  decreased 
when  ft  is  increased. 

Wait  and  Spies  [1964]  have  considered  the  case  of 
an  ionosphere  in  which  the  electron  density  increases 
upwards  in  an  exponential  model.  They  find,  in  this 
case,  that  the  effect  of  the  magnetic  field  is  to  cause  a 
reduction  in  the  phase  shift  on  reflection  for  propaga¬ 
tion  to  the  west.  Clearly,  for  the  admittedly  crude 
model  used  here,  this  should  produce  a  decrease  in 
in  D  for  propagation  to  the  west.  (This  is  confirmed 
by  the  second-order  mode  calculation  for  an  exponen¬ 
tial  ionosphere  in  a  supplement  to  Wait  and  Spies 
[1964].)  The  difference  in  the  effect  of  the  earth’s 
magnetic  field  on  the  reflection  phase  shift  for  the 
two  models  confirms  the  views  expressed  by  Wait 
and  Spies  [1964]  that  at  night  the  ionization  gradient 
is  quite  sharp.  Thus  the  behavior  is  perhaps  more 
like  that  of  a  sharp  rather  than  an  exponential 
boundary. 

Further  evidence  concerning  the  dependence  of  the 
phase  velocity  on  direction  of  propagation  is  contained 
in  a  paper  by  Swanson,  et  al.  [1965].  in  which  it  is 
shown  that  at  10.2  kc/s  the  diurnal  variation  on  the 
Hawaii-Cana!  Zone  path  is  about  10  percent  greater 
for  east-to-west  propagation  than  for  the  opposite 
direction  of  propagation.  Swanson,  et  al.  [1965J 
indicate  that  the  difference  occurs  at  night.  Thus 
the  nighttime  phase  velocity  for  east-to-west  propaga¬ 
tion  is  less  than  for  propagation  in  the  opposite  direc¬ 
tion.  Provided  that  this  difference  is  due  to  a  change 
in  phase  shift  on  reflection,  (2)  and  (3)  show  that  the 
phase  change  on  reflection,  ft,  must  be  greater  for 
east-to-west  propagation  than  fur  propagation  in  the 
opposite  direction.  This  is  also  in  agreement  with  the 
conclusions  reached  above. 


7.  Conclusions 

The  explanation  of  sunrise  and  sunset  fading  ob¬ 
served  on  long  VLF  paths  given  in  an  earlier  paper  is 
confirmed  by  data  for  propagation  in  the  opposite 
direction  to  that  initially  used.  Other  data  for  differ¬ 
ing  frequencies  and  propagation  direction  are  also 
considered.  It  is  found  that  the  difference  in  the 
phase  velocities  of  the  first-  and  second-order  modes 
under  nighttime  conditions  is  less  for  propagation 
towards  the  west  than  for  propagation  towards  the 
east.  This  is  found  to  he  qualitatively  consistent  with 
the  results  expected  for  a  sharply  hounded  ionosphere 
and  a  horizontal  terrestrial  magnetic  field. 
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Observations  on  signals  received  in  darkness  from 
the  east  show  that  it  is  not  necessary  mat  the  second- 
order  mode,  which  plays  an  important  part  in  the 
fading,  be  excited  by  the  transmitting  antenna;  the 
evidence  is  that  it  can  be  adequately  excited  at  or  near 
the  sunrise  terminator. 
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